Phasefield theory 
for fractional diffusion-reaction equations 
and applications 

Cyril Imbert*and Panagiotis E. Souganidis^ 
July 31, 2009 



Abstract. This paper is concerned with diffusion-reaction equations where 
the classical diffusion term, such as the Laplacian operator, is replaced with a 
singular integral term, such as the fractional Laplacian operator. As far as the 
reaction term is concerned, we consider bistable non-linearities. After properly 
rescaling (in time and space) these integro-differential evolution equations, we 
show that the limits of their solutions as the scaling parameter goes to zero 
exhibit interfaces moving by anisotropic mean curvature. The singularity and 
the unbounded support of the potential at stake are both the novelty and the 
challenging difficulty of this work. 
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1 Introduction 

This paper is concerned with diffusion-reaction equations where the classical 
diffusion is replaced with a singular integral term. Our aim is somewhat clas- 
sical: to show that the limit of their solutions after properly rescaling them in 
time and space exhibit a moving interface. However, we will deal with integral 
term whose potential is anisotropic, singular and with unbounded support. 
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9t^" + i - ^1 + fiun}=0 (2) 



Fractional diffusion-reaction equations. We consider for t > and x G 

R^, N>2, 

dtu' + -^!^-r^u' + f{u^)^ = o (1) 

and 

I In e 

where is a singular integral operator depending on a parameter a € (0, 2) 
and / is a bistable non-linearity. 

Example 1 (Standing example). The standing example for f andX^ are 

f[u) — u{u'^ — 1) and I^u — -^A^u 

where A "2 denotes the fractional Laplacian of the function u. We recall that 

dz 



A^u— {u{x + z) — u{x))- 



■A 



where the singular integral must be undertood in the sense of Cauchy's principal 
value. 

More generally, we will consider singular integral operators of the following 
form 

T^u{x) = J ^{x + ez) — u{x) — eDu(x) ■ z1b{£z)^ J{z)dz (3) 

where B denotes the unit ball and where the function J : R, which will 

be often referred to as the potential, can be of two types 

either J(^z) ^ g -1— or J e LI(R^) n C,(R^) (4) 

with a e [1,2) and g : S^^^ = {z e : \z\ = 1} ^ (0, +00) continuous and 
where Cc(IR.^) denotes the space of continuous functions with bounded support. 
The first potential will be referred to as the singular one while the second one 
will be referred to the regular one. As far as the standing example is concerned, 
J is singular with g = 1. 

Phasefield theory for diffusion-reaction equations. In [8], Chen proved 
rigourously that the solution of the AUen-Cahn equation yy generates a front 
moving with mean curvature as long as the front is regular. Thanks to definition 
of fronts past singularities [THl [5], Evans, Soner and the second author [T7] 
proved that this is still true after the appearence of singularities. Such results 
are generalized to a large class of bistable non-linearities by Barles, Soner and 
the second author [3] where a general phasefield theory for reaction-diffusion 
equations is introduced. In [T, T, an abstract method is developed in order 



2 



to deal with more general reaction-diffusion equations and to handle boundary 
conditions. In particular, non-local reaction-diffusion equations are considered 
in [?j but integral operators are not singular. As the proofs of the present paper 
will show it, it is a challenging difficulty to be overcome. 

Motivations. Recently, Caffarelli and the second author studied threshold 
dynamics-type algorithms corresponding to the fractional Laplace operator for 
a € (0,2). They proved that after properly rescaling them, they converge to 
an interface moving by mean curvature in the case a > 1 and to a fractional 
mean curvature in the case a < 1. Hypersurfaces with zero integral curvature 
are studied in ]7] . See also where the level-set approach [THl [5] is developed 
for such a geometric flow. 

As far as applications are concerned, two main physical models motivate 
the present study. The first application we have in mind is dislocation dynam- 
ics. Dislocation theory aims at explaining the plastic behaviour of materials 
by the motion of linear defects in crystals. Peirls-Nabarro models [55] consist 
in approximating the geometric motion of these defects by non-local diffusion- 
reaction equations such as ((2) . In [19] , such an approximation is also used and 
formal expansions are performed. In [20l |2T] , Garroni and Miiller study a vari- 
ational model for dislocations that can be viewed as the variational formulation 
of the stationary version of ^ . 

The second application we have in mind is statistical mechanics and more 
precisely stochastic Ising models. These models were introduced by Kac, Uhlen- 
beck and Hemmer p4] (see also [29]) to justify the validity of the Van der Waal's 
phase diagram. The interaction beween particles is described by the Kac poten- 
tial. A lot of work has been done since then to understand the hydrodynamic 
limits of such interacting particle systems and it is beyond the scope of this 
paper to give a complete list of references. However, we can mention the papers 
by De Masi, Orlandi, Presutti and Triolo [IS [H [H [H [H] and Katsoulakis 
and the second author [25l [26l [27] . The interested reader is also referred to the 
monograph of De Masi and Presutti [TB] and the book of Spohn [3D|. In the 
papers we mentioned before, hydrodynamic limit of stochastic Ising models with 
general dynamics are studied. In particular, a mean field equation is derived, 
see |12| . By many ways, these equations can be viewed as non-local reaction- 
diffusion equations. The next step is to show that for appropriate scalings the 
solution of the mean field equation approximates an anisotropic mean curvature 
motion; see for instance [111 I25j . Green-Kubo type formulae are provided for 
the mobility and the diffusion matrix in terms of a standing wave associated 
with the mean field equation. 

On one hand, the Kac potential is assumed, in most papers, to be reg- 
ular with a compact support. On the other hand, Lebowitz and Penrose |29i 
Eq. (1.20b), (1.21a), p. 100] consider potentials J that are singular; more precisely, 
they assume that for small z, the singularity of J is of the form for 
a > 0. 
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Description of the results. Our main result states that, as e ^ 0, the so- 
lutions of ll]) and ([2]) can only have two limits: the stable equilibria of the 
bistable non-linearity / (see Section [5] for definitions). The resulting interface 
evolves by anisotropic mean curvature; moreover, Green-Kubo-type formulae are 
obtained: the mobility and the diffusion matrix of the geometric flow are ex- 
pressed in terms of the standing wave associated with the bi-stable non-linearity; 
see Eq. ((22)) . ((23)) and (|25|) below. Even if the proof follows the classical idea 
of constructing barriers by using traveling waves, the reader will see that classi- 
cal arguments fail when extending the barrier away from the front; several new 
ideas are needed to handle the unboundedness of the support. In order to handle 
the anisotropy of the potential, we have to use ideas developed by Katsoulakis 
and the second author [2^ and introduce correctors to cancel oscillating terms 
by averaging them. This implies in particular that anisotropic traveling waves 
must be considered. But because the integral term involves a singular poten- 
tial, passing to the limit in averaged oscillating terms is challenging and this 
constitutes the core of the proof of the convergence theorem. 

As the reader will see it when going through the preliminary section or in the 
statement of the convergence theorem, several assumptions on traveling waves 
and the linearized traveling wave equation are necessary (if not mandatory). 
Even if we do not construct such waves and correctors and assume that they 
exist, the reader can check that the assumptions we make are natural. For 
instance, the decay estimate (|13p is expected since its corresponds to the one 
of the kernel of the fractional Laplacian in the one dimensional space. See also 
[6]. We plan to construct them in a compagnion paper. 

As explained above, we will consider two kinds of potentials: singular and 
regular ones. As far as the singular case is concerned, we distinguish two sub- 
cases, depending how singular is the potential at the origin. Since potentials 
are positively homogeneous in the singular case, potentials in the subcase a = 1 
decay as when \z\ — > -|-oo. This corresponds to the dislocation dynam- 

ics model. As the reader can see it, the scaling involves a logarithmic term; this 
factor is well-known in physics and the interested reader is referred to [5] for 
instance; see also [TUJIB]. An additional comment about singular and regular 
potentials concern the Green-Kubo-type formulae. It turns out that these for- 
mulae are different in singular and regular cases. However, we give in appendix 
a formal argument to shed some light on the link between these two formulae. 



Additional comments. As the reader can see it, we are not able to deal 
with the case a < 1 even if, in view of the results of [6], we should observe an 
interface moving with fractional mean curvature (see Section [2] for a definition). 
In this case, the equation should be rescaled in time as follows 

dtu^ + ^^-I^u^ + f{u^)^ (5) 

with 



X^u{x) ~ / [ u{x + ez) — u{x) \ g 



z \ dz 
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for some a € (0, 1). The reader can check that we are able to pass to the Umit 
in (the average of) oscillating terms (see Lemma [10] below), which is usually the 
difficult part of the convergence proof. We are even able to construct a barrier 
close to the front. But because the diffusion-reaction is non-local, we are stuck 
with extending the solution away from it. In particular, the very slow decay of 
the potential at infinity does not permit us to use the new ideas we introduced 
in the singular case a > 1. This difficulty is unexpected since in [6], this case is 
the easiest one. We hope to find a path toward this result in a future work. 

In the one dimensional space, moving interfaces are points. Gonzalez and 
Monneau [52] considered such a case and proved a result analogous to our main 
one by taking advantage of the fact that the limit is a (system of) ordinary 
differential equation(s). In particular, the restriction on the strength of the 
singularity can be relaxed in this case. 

Organization of the article. The first section is devoted to preliminaries. 
In particular, traveling waves are introduced as well as the linearized traveling 
wave equation which is the equation satisfied by correctors; see Subsections 12.31 
and 12.41 We also introduce the geometric motion by mean curvature (Subsec- 
tion 12. 5p together with its equivalent definition in terms of generalized flows 
(Subsection 12. 6p . Our main result is stated in Section [31 In the remaining of 
this section, we explain how to reduce the proof of this convergence result to 
the construction of an appropriate barrier (see above). Section [3| is dedicated 
to this construction. The last section (Section [5]) contains to core of the proof 
of the convergence result: the limit of the average of oscillating terms. Finally, 
we give in appendix a formal argument to explain the link between the two 
Green-Kubo formulae obtained in the convergence theorem. 

Notation. The Euclidian norm of a; e is denoted by The ball of 
center x and of radius r is denoted by B^(x). We simply write for -8^(0) 
and B = Bi denotes the unit ball. The scalar product of x and y is denoted by 
X ■ y. The unit sphere of K.^ is denoted by §^~^. The set of symmetric N x N 
matrices is denoted by Sn- The identity matrix (in any dimension) is denoted 
by/. 

Given two real numbers a, 6 g R, a V 6 denotes max(a, b) and a Ab denotes 
min(a, b), denotes a V and denotes —(a A 0). In particular, > and 
a = a"*" — a~. 

Given a set A, 1a denotes its indicator function that equals 1 in A and 
outside. The signed distance function d^ix) associated with A equals the 
distance function to R^ \ A if x e A and the opposite of the distance function 
to A if X ^ A. 

The set of continuous functions / : R with compact support is denoted 

by C°. 

Given a family of locally bounded functions : C R'' ^ R indexed by 
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e > 0, the relaxed upper and lower limits are defined as follows 

liniinf ^/^(x) = liminf f^{y) and limsup ^/^(a;) = limsup . 

fe = f for any e > 0, these relaxed semi-limits coincide with the lower and 
upper semi-continuous of a locally bounded function /. 

For traveling waves q{r^ e) and correctors Q{r, e), q and Q denote derivatives 
with respect to r. 



2 Preliminaries 

This section is devoted to the presentation of the assumptions we make about 
non-linearities, traveling waves and linearized traveling wave equations. We 
also briefly describe the construction of fronts at stake after rescaling fractional 
diffusion-reaction equations. 

2.1 Fractional diffusion-reaction equations 

We can write ((1]), ^ and (O as follows 

dtu' + j^(^~r^u' + f{u')'^=0 (6) 

with 

{£ if a > 1 , 

e| Inel if a = 1 , (7) 
e" if a < 1 

(see Remark H}. 

We will use later on that the potential J satisfies in the singular case the 
following properties 

' J is smooth on \ {0}, even and non-negative 
< \J{z)\<j^ (8) 
Ji^) ^ 9 (]f|) as \z\ +00 

with a e (0, 2). We also mention that if a < 1, then 

\z\'Mz) + \V{\z\'Mz))\ < K{z) e L\B) (9) 
where J,{z) = e-^-°'J{e-^z). 

2.2 Bistable non-linearity 

We make the following assumptions. 
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Assumption 1 (Bistable non-linearity). The non-linearity f : M. ^ W is 

and such that 

• for all h G (0,H), they are constants m±{h) and mQ^h) such that 

f{mi{h))=h, i£{0,+,—} and m-{h) < mo{h) < m^{h). 

(10) 

• / > m (m_,mo) and f < in (rfiQ^rh^). 

• f'{rh^) > and /'(mo) < where fh± = m±{0) and fho — mo(0). 
In particular, if h > 0, we have fa± < m±{h) and fho > mo{h) and 

\m±{h)-ffi±\<Cfh. (11) 

2.3 Anisotropic traveling wave 

In this subsection, we describe the anisotropic traveling waves we will use in 

the construction of barriers in order to get the main convergence result. In par- 
ticular, we make precise the decay we expect for such waves. This construction 
will be achieved in a future work. 

Assumption 2 (Anisotropic traveling wave). For h £ {0,H), there then exist 
two continuous functions q :Rx S^~-^ — > M and c : M x S^~^ R such that 

q{r, e, h) — > m±{h) as r ^ ±oo 

(the limit being uniform with respect to e G S^~^) and 

cq-Ie[q] + f{q)=h (12) 

where 



/ + e-z)- q{0 - qiOe ■ ^1b(2) j J{z)dz if a > 1 
I [q{^ + e-z)- 9(0) J{z)dz ifa<l 

for any e G E>^~^ . 

The traveling wave q is increasing in r and the following estimates hold true 
for any (r,e) e M x S^-^ 

\q{r,e,h)-m±{h)\ = O (^j-^^ as r ^ ±oo (13) 

SUp{||De«||oo + Pe%9||oo} < (14) 
h>0 

for some constant C > and with a limit uniform in e G S'^~^. The function q 
also satisfies 
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as h and the limit is uniform in e. 
The speed c satisfies 

hc{e, h)>0 if h^O (16) 

and 

^^^c(e) as h^O (17) 

and the previous limit is uniform with respect to e. Moreover, the function c(e) 
is continuous on S^^^. 



Standing wave. Notice that ([TT]) implies in particular that c(e, 0) = 0. Hence 
q{e, 0) is a standing wave. It is denoted by g*^ in the remaing of the paper. 



Reduced integral operator. The operator does not depend on e if J{z) 
is radially symmetric, i.e. when J{z) = j(|2;|). We illustrate this fact in the 
next lemma where Te is computed in the case where J{z) = g(i)|z|~^~". 

Lemma 1. Assume that J{z) — g(pj) \z\n+c. ■ Then 



r) = an(e) / + r) - g(0)i„i+„ 



where 

I \ f du 



(1+ |u|2)(W+")/2 



Remark 1. We recognize the fractional Laplacian of order a in the one dimen- 
sional space (up to a multiplicative constant). 



2.4 Linearized traveling wave equation 

In this subsection, the linearized traveling wave equation is considered. 
Loosely speaking, we need to know that the kernel of the linearized oper- 
ator C reduces to Mg and, if / is regular enough, so is the solution Q of 
^QiO — Pfi^i^i^) where Pf is the projection of / on the space orthogo- 
nal to q. We need in particular to be able to say that Q decays at infinity. Let 
us be more precise now. 

The linearized operator C associated with ()12p around a solution q is 

CQ^cQ-I,[Q] + f {q)Q. (18) 

Given a smooth function d : [0, T] x M, let us consider 

I , ^ ^ f f t ^ dit,x + ez)-dit,x) At( \^ 

a{r,e,t,x) = J / < g(r H ,e) — q[r + e ■ z,e) >J(z)dz 

where h = e (resp. e\ lne|, e") if a > 1 (resp. a = 1, a < 1). 
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Assumption 3 (The linearized TW equation). 

Ker L = Ker (£)* = span q. 

Moreover, if d : [0, T] x ^ R is a smooth function, there then exists a 
continuous solution Q:Rx S^-^ x [0,r] x to 

CQ = Pga{-,e,t,x) 

where Pq stands for the projection on the space orthogonal to span q. In partic- 
ular, there exists Cq > such that for any h, e, ^, 

\Q\ + \dtQ\ + \D,Q\ + \D,Q\ + \Dl,Q\ + \Dl,Q\ < Cq (19) 

Q^O asr ^ +oo (20) 

where Cq does not depend on h, e, t, x and the limit is uniform in h, e, t, x. 
2.5 Geometric motions 

In this subsection, we introduce the geometric motions of fronts at stake 
when rescaUng the fractional diffusion-reaction equations. 

It is well-known that singularities can appear on the front in finite time 
when considering, for instance, the mean curvature motion. We thus classically 
use the level-set approach to define a front for all times. We recall that this 
approach consists in looking for a front Ft under the form {x : u{t, x) = 0} and 
to derive a PDE satisfied by u. 



Anisotropic mean curvature motion. In the case of an anisotropic mean 
curvature motion, we obtain the following degenerate and singular parabolic 
equation, 

dtu = n [pv^ Tr ((J - O Du)A [Dv^ D'^u^ (21) 

where /x : §^~^ — > M+ and A : S^~^ 5jv are continuous functions and I 
stands for the N x N identity matrix and e = . We will see that the function 
fj, (which will be referred to as the mobibility) is given by the following formula 

Me) = {y'(9")'(e,e)de| (22) 

As far as the function A is concerned, we distinguish cases. In the singular case 

and if a > 1, we have for all e e E>^~^ 

^(e) = |/ l{q°{^ + z^,e)-q°{^,e)rj^d^y,{e) (23) 

with 

A,{e) = a{a - 1) J^^_^ (1, «) ® (1, u)9{l, ^ |^|2)(^+»)/2 (24) 
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(where the space orthogonal to e is identified with ^). If a = 1, 

(m+-TO_)2/ (9(»9J{e)]a{de) . (25) 

In the regular case, we have for all e G S^^-"^ 

A(e) = I qO{^,e)q°{£,+e- z)z(»zJ{z)dz. (26) 



Remark 2. As a matter of fact, in the singular case with a > 1, ^'s given by 
and (|26p are the same, at least formally. But it is not even clear that 
the integral defining A{e) in is well defined in the case a > 1. A formal 
argument is given in Appendix. 

Fractional mean curvature motion. A fractional version of this motion 
can be defined. More precisely, for any a < 1, one can consider the following 
PDE, 



dtu = fi[Du] K[x,u{t,-)]\Du\ (27) 
where fj, is defined by and 



k[x,U] ^ K*[x,U] ^ I \l{U{x+z)>U(t,x),e-z<0} 

z \ dz 



\N+a 



— '^{U(t,x+z)<U{t,x),e-z>0} \9 

This can also be written under the general form 

K*[x,U] = v{z e : U{x + z) > U{t,x),e-z < 0} 

- i^{z e : U{t, x + z) < U{t, x),e- z>0} 

for some non-negative Borel measure v which is eventually singular. A general 
theory is developed in 23J to prove that the geometric fiow is well defined. The 
definition of a viscosity solution for ([?7|) implies the use of the following quantity 

^ \^x , U] = v{z £ : U{x + z)> U{t, x), e • z < 0} 

- v{z £ : U{t, x + z)< U(t, x),e- z>0}. 

Geometric non-linearities. In the following, we will use the notation: if 
a > 1, 

Fip,X) = -A/(p)Tr (i(p)A) - c(0,j5)b| (28) 

and if a < 1, 

F(p, X, [</.] ) = -fiip)K* [x,^]\p\- c(0, p) \p\ (29) 
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for p ^ and p — p/\p\. Since non-linearities F are discontinuous, it is necessary 
to use the lower and upper semi-continuous envelopes and F* of F in order 
to define viscosity solutions of (|2ip . In the case a < 1, we have 

F*{p,x, [4>]) = -^lMK,[x,4>]\p\-UQ,p)\p\ 

F, {p, X, [cj)] ) = -fi* {p)k* [x,^]\p\~^ c* (0, p) \p\ 
with the convention — 0. 

2.6 Generalized flows 

As we explained it above, the level-set approach is necessary in order to 
define the anisotropic mean curvature motion of a curvature after the onset 
of singularities. It is proved in [4J (see also [2]) that this notion of solution 
is intimately related with the notion of generalized flows of interfaces whose 
definition is recalled next. 

Definition 1 (Generalized flow for pT|) ). An family — {ilt)t>o (resp. U = 
{Ut)t>o) of open (resp. closed) sets o/ is a generalized super-flow (resp. 
generalized sub-flowj of (j21|) if for all {to,xo) G (0, +oo) x M^, h > 0, and for 
all smooth function cj) : (0; +00) x — > R such that 

(i) (Boundedness) There exists r > such that 

{{t, x) e [to, to + /i] X : (j){t, a;) > 0} C [to, + h] x B{xo,r), 

(ii) (Speed) There exists > such that 

dt<j) + F*{D(l),D'^(j)) < in [tQ,tQ + h] x B{xo,r) 
(resp. dt(l) + F^{D(t), D^cj)) > in [to, to + h] x B{xo,r)), 

(Hi) (N on- degeneracy) 

D(t)^Q in {(s, y) e [to, to + h]x B{xo,r) : (j){s, y) = 0}, 

(iv) (Initial condition) 

{yeB{xo,r) : (f>{to , y) > 0} C , 
(resp. {y G B{xo,r) : (j){to,y) < 0} C \ Ft J, 

then 

{y e B{xo, r) : 0(to -f /i, y) > 0} C nt,+h 
(resp. {y G B{xo, r) : 0(to + /i, y) < 0} C R^ \ F^.+n)- 

Remark 3. Remark that this definition slightly differs from the one introduced 
in [4]. However, a quick look at the proof of the abstract method from [4J that 
will be used below should convince the reader that this definition is adequate 
too. 
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3 The convergence result 



This section is devoted to statement of the main resuh of this paper. We also 
explain why its proof reduces to the construction of an appropriate "barrier" 
which will be constructed in the next section. 

Loosely speaking, we will prove that solutions of the fractional diffusion- 
reaction equation ^ approximate, as e go to 0, the motion of a front moving 
with a normal speed equal to its mean curvature. Moreover, the results state 
that the mean curvature motion is anisotropic and that mobilities and diffusion 
matrices are given by Green-Kubo formulae (see (P^ . and (HH)). 



3.1 Statement of the main result 

We now state our main result. 

Theorem 1 (Convergence result when a > 1). Let J be given by (|4]) with 
a > 1 in the singular case and let f be a bistable non-linearity. We suppose 
that Assumptions]^ [H are satisfied by f. Let u"^ be the unique solution of 
^ if a > 1 and of ^ if a = 1 associated with a continuous initial datum 
Uq : [to_, ?7i-|_] defined by 

ut,{x)^q'>(^,Ddo{x)) (30) 



where is the standing wave associated with the diffusion-reaction equation 
and do is the signed distance function to the boundary of a smooth set Qq. 

Let u be the unique solution of the geometric equation pil) supplemented 
with the initial condition u{0,x) — doi^x), where fi is given by ()22|) and A is 
defined in (j26p in the regular case and (|23p . (|25p in the singular case. 

Then the function satisfies 

— !■ TOj- in |u > 0} 
£ - • r ^ ni as e 

where fh± denote the stable zeros of f ; moreover both limits are local uniform. 

Proof. In order to prove this theorem, we use the abstract method developed in 
[4] and [2]. Consider two open sets defined as 

= Int{(t,.x) e (0,+oo) X : liminf. " ~ "^^ > 0} C (0, +oo) x 

V 

= Int{ (t, x) e (0, +oo) X R" : lim sup * " ~ < 0} C (0, +oo) x R^ 

V 

where interior is considered with respect to (0, +00) x R^. 

We next define their traces at initial time by considering the lower semi- 
continuous function x = lf2i ~ and the upper semi-continuous function 
X = l(n2)c — lf22 defined on (0, -l-oo) x R^. They can be extended at i = 
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by setting x(0,a:) = liminft^o.y^x y) and x(0,a;) = limsupj^g 2/)- 
We now define 

nl^{xeR^ : x(0, x) = 1} and nl = {xem.^ : x(0, x) = -1} . 

The metliod developed in [4j| consists in proving the following propositions. 

Proposition 1 (Initial time). The set {x e : do{x) > 0} is contained in 
JIq. Similarly, the set {x G : dQ{x) < 0} is contained in 57q. 

Proposition 2 (Propagation). The set il^ (resp. ^'^) defined above is a gen- 
eralized super- flow (resp. sub- flow) of (|2ip . 

The proofs of both propositions are postponed. Applying next [21 Corol- 
lary 2.1] (see also |4] Corollary 3.1]), we conclude the proof of Theorem [TJ □ 

It remains to prove Propositions [1] and [21 Both rely on the construction of 
barriers for smooth fronts; in our case, the term "barrier" refers to a sub- or 
super-solution of the fractional diffusion-reaction equation. 

3.2 Proofs of Propositions [H and [2l 

As we shall see it, proofs of Propositions [H and [2l reduce to the proof of the 
following one. 

Proposition 3 (Construction of a barrier). Given to > and xq G M.^ , con- 
sider a smooth function cj) : (0, -|-c») x M.^ — > R such that (i) , (ii) (Hi) from 



Definitions^ are satisfied. For all P >0, there exists a sub-solution If^'^ of ([T]) 
if a > I and ([2]) if a — 1 such that 

U'^^ih^x) < (™+(-/3r/) -/3r;)l{d(,„..)>^j +m_(-/377)l{d(t„..)<^} for x G 

(31) 

where d{t, x) denotes the signed distance to the set {y : (j}{s, y) — 0} which has 
the same signs as <j) and m±{—Pri) are the stable equilibria of f + /3ri. Moreover, 
ifd{t,x) > —2/3, then 

liminf , -{t,x) > -{Cf+2)f3 (32) 

where Cf appears in (jlip . 

We first derive Proposition [21 from the construction of the barrier. 

Proof of Proposition [H Let us explain why Proposition[31together with the com- 
parison principle for ([6]) yield the desired result. 

We consider a smooth function such that (i)|(iv) hold true. Let d denote 



the signed distance function to {0 = 0}. We derive from (iv) that 



{dito, •) > 0} C {0(io, •) > 0} C nl = {hminf ^ > 0} 
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For the sake of clarity, m± denotes rn±(~ (irj) . And we will do so in the remaining 
of the paper. Hence 

— Til ^ 

{d{to, •) > /3} C int{liminf , — > 0} . 

1] 



Hence by using (i) we conclude that it^ satisfies on one hand 

> m+ — /3?7 > m+ — Pt] in {d{to, ■) > 13} . 

On the other hand, since to_ is a trivial solution of the diffusion-reaction equa- 
tion, we have > rh- > m„. We thus conclude that satisfies 

w^(to, x) > {m+ - /5?7)l{d(to,-)>/3} + ™-l{rf(to,-)</3} for x e , (33) 

We now use Proposition [3] in order to get a sub-solution with the desired 
properties. Combining ((3T|) and ((33)) yields that If^ < u"^ at t = to. We thus 
conclude by using the comparison principle for ([6]) that < u'^ on [to, to + ft.] x 
and ^ implies that 

liminf. "^(^° + ^^"^''^+ > -(Q 2)/3 

as soon as (i(to + h,x) > 2(3. Since /3 is arbitrary, the proof is complete. □ 

We now prove Proposition [TJ 

Proof of Proposition QJ We only prove the result for fig since the proof for 
is similar. 

Let xo be such that dQ{xo) =: 25 > 0. We have to prove that xq G fip. In 
other words, for all {t,x) in a neighbourhood of (0, xo), we would like to prove 

lim mi * ( t, x)>0. 

V 

In order to get such a result, we construct for any small /5 > a subsolution 
U^^P of ® such that 

U^'^{0,x) < u%0,x) 

and satisfying ([32]) for some function d{t, x) such that {d > 2(3} contains a 
neighbourhood of (0,a;o). 

There exists r > such that for any x £ B{xo, r), do{x) > S > 0. Consider 
next the smooth function 

(bit,x) - {r - Ct)l - \x - xo\\ 

The associated distance function is given by the following formula 

d{t, x) = r — Ct — |a; — a;o| . 

Remark that {d > 2(3} = Ut>o{t} x B{xo, r ~ Ct - 2(3). 
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We claim that ([55)1 holds true with to = 0. Indeed, when d(0, x) > (3, we 
know that ^0(2^) > 5 and this implies that u^(0, x) > m^ — ^rj for e small enough 
as showed now 



u^{Q,x) 



> 



do{x) 



,Dda{x] 



> q"[-,Ddo{x] 



> 
> 
> 



+ o(ei+") 
+ o(ei+") 



Notice that (i) and (iii) are satisfied. As far as |(ii)| is concerned, it is only 
used in the construction of the barrier in order to get (|34p below. Wc thus have 
to prove that we can choose C > such that ([M)) also holds true. The constant 
C is chosen as follows 



C > sup (-/i(e)Tr(A(e)) - c(e)) 



and ((34)) holds true for 7 and h small enough. 



□ 



4 Construction of the barrier 

This section is devoted to the proof of Proposition [S] 

Proof of Proposition[3i The proof proceeds in several steps. We first construct 
a sub-solution U of the diffusion-reaction equation close to the smooth front and 
we then extend it to the whole space. 



A barrier close to the front. Using (i) (ii) and (iii) we know that there 



exists 7 > such that d is smooth on the set 

Q7 = {Ml < 7} 

and D(j){s^x) ^ on and 

dtd < fi{Dd)TiiA{Dd)D^d) - inQ^. (34) 

We used the fact that \Dd\ = 1 in which also implies that D^dDd = in 
Q~f. For (3 < 7/2, we next define a "barrier" as 

U{t, x) = q( ^^h^l^^,V{t, x)-(3h 



hQ I — '^,V{t,x),t,x,-(3h] -2(3h 
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where h > will be chosen later, q denotes the traveling wave given by As- 
sumption [21 the function V is assumed to be smooth, to coincide with Dd in 
and to be such that ^ < l^?! < |. Let us point out that we would like to 
choose V = Dd but this function is not well defined everywhere away from the 
front and even if we prove that C/ is a subsolution close to the front, U has to 
be defined everywhere since the diffusion-reaction equation is not local. As far 
as the function Q is concerned, it will be chosen later. 

Plugging the barrier into the diffusion-reaction equation. In order to 
prove that the barrier we introduced in the previous step is a sub-solution of the 
diffusion-reaction equation close to the front, we first plug it into the equation. 

Lemma 2. If h — j] and [3 < [3 (depending only on 5^), then the function U 
satisfies the following inequality in 



dtU + -{~TJJ + f{U)} < 

er] e 



M(e)Tr(^(e)B) - f^,a, - ^ 



e 



+ - [£Q{r,e,t,x) - {ae[r,e,t,x) - q{r,e)fisas)] + (err) (35) 

£ 

where r — ^^ih^lzlS.^ g — Dd{t,x), B = D'^d{t,x), ^ctO^c & R two real 
numbers to be chosen later and 

a,{r,e,t,x)^^J |g(r+ + ~ , e)-qir + e-z, e)| J(z)dz . (36) 

As far as error terms are concerned, we have 

(err) = ^R[T,] + -R[Tq] + - f-g^ - 2f{q)f3) + o{e-') (37) 
£ry £ £\8 / 

and 

^^^'^^^^^^'^^ ,V{t,x + ez) ] +BDeqir,e) ■ ezleiez) j>J(z)dz, (38) 



R[Tq] = J I^Q{r + e- z,e,t,x) 



^ 'd(t,x + ez)-2l3 ^, 
Q ( -,V{t,x + ez),t,x + ez 



{BDeQ + D^Q) ■ ezlgiez) \j{z)dz . (39) 
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Proof. We compute the quantity dtU + (stj)-'^ {-I^U + fiU)}. 

The fact that g is a traveUng wave (see (|T2|) ) together with a uniform bound 
on Q with respect to all its variables (we will choose Q below so that it satisfies 
such a condition) permits to get 

f{U) = f{q) + hf'{q)Q-2f'{q)(3h + 0{h') 

= -0h +J,[q] - c{e,-ph)q + hf'{q)Q - 2f'{q)(3h + 0{h'') . 



Rearranging terms, we thus obtain, for (t, x) e 



c(e, —Ph) 
V 



{0,Dd) 



ei] erj 



+ -(Xe[<7] + /i/'(g)Q-ISC/) + 0(-) (40) 
er] erj 



We immediately see from this computation and in view of (jl7p that h must be 
chosen as follows 

h = Tj . 

We next write 



XeM(r, e) + ilf'{q)Q ~ IlU{t, x) = rjTg + T, 



where 



and 



Tq = /' (g) Q - lUQie-'id - 213), Dd, t, •)) = CQ - cQ + R[Tq] 

T, - l4q](r)~I'Me-Hd-2f3),V)]{x) 

q{r + e ■ z,e) — q{r, e) — g(r, e)e • zlB{z)^J{z)dz 
( d{t,x + ez) - 2(3 



(41) 

(42) 



-,V{t, X + ez) ] ~ q{r, e) 



(^q{r, e)- + BD^q{r, e)j ■ ezlsiez) >J{z)dz . 



Hence, T„ can be written as follows 



Tq = -rja^ + R[Tq 



(43) 



We now compute the time derivative of the barrier U . We use (j34p in order 
to get 

dtU = i [dtd] + {D,q + jjDeQ) ■ D{dtd) + rje-^Qdtd + r^dtQ 



<1 



j,{e)TriA{e)B) ~ f 



{D,q + rjD.Q) ■ D{dtd) + Tje-^Qdtd + 7]dtQ. 



(44) 
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We next combine (gOl), dH]), gH), (gS]) and dH]) to get 



9t{/+^{-XSf/ + /(£/)}<! 



n{e)TT{A{e)B) - — — -- 

+ -[-a, + £Q] + (err) 
e 



with 



(err) = i-i?[r,] + Ji?[rQ] - + ^ (-^^ - V{q)f3 

+ {D,q + ijD.Q) ■ D{dtd) + i^e-^Qdtd + rjdtQ + o{e-^) . 
Using P?]) . and p^ . we finally get ([55]) with the associated error term. □ 

Estimating error terms. In this paragraph, we prove that the right hand 
side of (j35p is non-positive. We first construct a corrector Q in order to handle 
oscillating terms. 

Lemma 3 (Choice of the corrector Q). There exist /ie,ae € M such that there 
exists Q satisfying 

CQ = flg — qUeO-e ■ (45) 
Proof. In view of Proposition [31 it is enough to choose /Zg and ds such that 

\a,{0-ti,d,q{0)q{Od^^O. 

The following choices permit to ensure such a condition 



and 

de{e, t,x) ^ J g(^, e)ae{^, e, t, x)d^ 

= ^ J J g(C, e)|g(^ + e • z + eW{t, x,z), e) - q{^ + e ■ z, e) )'J{z)d£,dz 

with 



VF(i, a;, 2:) = \d(t, x + ez) — d(t, x) — eDd(t, x) ■ z] 



□ 



Remark 4. The choice of h when rescaling fractional diffusion-reaction equations 
([6]) is made such that dg has a limit as £ ^ 0. 

The following lemma is the core of the proof of Theorem [T] and its proof is 
rather involved. This is the reason why we postpone it until Section [S] 
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Lemma 4 (Uniform convergence of approximate coefficients (I)). As e — > 0, 

deie,t,x) -> tr {A{e)D'^d{t,x)) 

and the limit is uniform with respect to {e,t,x) £ x Q^. 

We next treat error terms appearing in (err). 
Lemma 5 (Error terms (err)). We have 

R[T,] = o(e") = 0(7?) (46) 

and 

R[Tq] = 0(1) (47) 
uniformly in (e, t, x) E S^^^ x and for all r E M. and (3 < (3 — (3{S^) 

-g(r)^-2/'(g(r))/3<0. (48) 
Proof. We first prove Through a change of variables, we get 

^q(^^^^^^-^^^—^,V{t,x + z) ] ~BD,q{r,e)-zlB{z) )J,{z)dz 



where Je{z) = £-(^+") J(e"iz). By using ([5]) and ITU) , dominated convergence 
theorem permits to conclude. 

We next turn to the proof of (|47|) . To prove it, we first write 



R[Tq]^R^[Tq\-R^[Tq] (49) 

with 

R\Tq] - y|Q(r + e.z,e,t,a;)-Q(^^^^l^i|^^^,e,i,a;) |j(z)dz, 
R^[Tq\ = J I^Q{r + e- z,V{t,x + ez),t,x + ez) 

-Q{r + e ■ z,e,t,x) ~ {BD^Q + D^Q) ■ ezlB{ez)^ J{z)dz . 
As far as R^ [Tq] is concerned, we can write for any R> 
R'[Tq] < ||Q||oo^||i?'d||L-(s(..si?))£ / |z|V(z)dz 
+2||Q||oo / J{z)dz 

J\z\>R 

< C\\Q\\ool\\D^d\\L^^Bi.^,R))£R'-'' + 2C\\Q\\^R-" 
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where we used ^ to get the second inequahty. Choose now R such that eR<l, 
R +00 and £i?^~" 0; for instance R = e^^l"^ permits to conchidc in this 
case. 

As far as i?^[TQ] is concerned, we use once again ([S]) in order to write 



|i?2[TQ]|<e" 

- g(. . . , V{t, x),t, x) ~ {BD,Q + D,Q) ■ z1b{z) 



(. . . X + z),t,x + z) 

dz 



C. 



\N+a 



where C depends on sup^{||(5||oo + \\D1_eQ\\oo + II^^.^QIIoo} that is bounded by 
assumption (see Estimate ([T9|) ). 

It remains to prove (j48|) . It is enough to prove that there exists a constant 
Ctw which does not depend on h and such that for all r G R 

q{r) + Ctn,f'{q{r))>Q. (50) 

This inequality is trivial when f'{q{r)) > 0. Hence, we consider r such that 
f'{q{f)) < 0, that is to say 

m_ < q- < q{r) < q+ < fn+ 

for some constants q± which do not depend on h. If r satisfies the previous 
inequality, we deduce from (|13p that 

\r\ < R 

for some constant R which does not depend either on h. Now ([50]) is clear. It is 
enough to find an estimate from below for q on [—R, R] which does not depend 
on h. The proof of the lemma is now complete. □ 

Using Lemmata [3l H] and [5] we derive from ([35| the following inequality 

d,U + -{-KU + fm < ^ + oi-) < . (51) 
er] e e 2e 



Extension of the barrier away from the front. The remaining of the con- 
struction of the barrier consists in extending the subsolution U we constructed 
before in order that it is a subsolution in [to, to + ft,] x B(xo,r) (in particular, 
far from the front). More precisely, we modify U far from the front. Following 
we proceed in two steps. We first extend it on {d < 7} by to_ and then 
extend it on {d > 7} by m+ — (3r] The difficulty is to keep it a subsolution. We 
do this by truncating properly U. Truncating it from below by m_ is easy but 
truncating it from above by — (3r] is more delicate. 
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Upper estimates for U. We start by estimating from above the "barrier" 
function U we constructed before. We claim that the foUowing inequahties hold 
true 

U{t,x) < m_ in{d</3}, (52) 
U{t, x) < m+ - f3r] . (53) 

We first justify ([5^ . In view of the definition of C/, we use and (^0]) . in 
order to get 

U{t, x) <m^+ 0{e^+") - 2(3ri < ru- . 

We next justify ([55)1 by adapting an argument from First, (|20p implies 
that there exists c > which does not depend on h such that we have for |r| > c 

\Q{r,e,t,x)\<(3. 

Next, we claim that there exists i'{c) > such that we have for \r\ < c 

q{r) < 171+ — v{c) . (54) 
Now if \dit, x) - 2/3| > ec, then 

U {t, x) < 171+ + I3r] — 2/377 ^ ~ ■ 
In the other case, \d{t,x) — 2(3\ < ec, then 

U{t, x) < m+ - iy{c) + ryllQIloo - 2/3?; < m+ - P?] 
as soon as ?7|lQ||oo < '^(c)- 

Definition of U. We define for {t,x) e [to,to + h] x 

U{t,x) = max(C/(t, x), m_) . 
From dnil) and we get 

t7 = m_ in {d < /3} , (55) 
m_ < U < m+ — Pi] . (56) 

On one hand, a classical argument implies that C/ is a subsolution of ([6]) on 
= {—7 < d < 7} since it is the maximum of two subsolutions. On the other 
hand, (|5^ implies that U{t,x) — m_ on {d < —7/3}. Thus (7 is a subsolution 
on {d < 7}. 

We also shed light on the fact that U satisfies ([51]) at points of where 
U = U. This will be used later on. We will also need the following lemma. 

Lemma 6 (Gradient estimate for the barrier). There exists C > such that 
for all {t, x) e [to, to + h] x M.^ and all e > small enough, 

\eDU{t,x)\ < C . 
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Proof. This is a simple consequence of the definition of U and of Estimates (I14p 
and (HHI. □ 



Definition of W^'^ . We finally define 

Ue,p^l ^(rf)C/+(l-^(d))(™+-/3r?) ifrf<7, 
[ m+ — /377 if (i > 7 

where t/; : R ^ [0, 1] is a smooth function such that ip{r) = 1 if r < 7/2, ip{r) = 
if r > 37/4. We will see below that it is convenient to assume additionally that 
ip{5'-f/8 + r) — 1 ~ ^(57/8 — r). We deduce from properties of tp and (ISS)) and 
that 

?7^''' = m_ in {d < P} (57) 

J/^''^ = m+ - /3?7 in {d > 87/4} (58) 

m_ < C/^^^ <m+- Pr]. (59) 



In particular, (|3T|) is clearly satisfied. 

We also deduce from ([55)1 and the definitions of U and C/^^'' that 

C/ < [7 < U'-'^ . (60) 

In particular, if d{t,x) > 2/?, then > 2/? in a neighbourhood V of (t, x). In 
particular, 

ljs,p >u = m+ + o{t]) ~ 2/3?7 in V . 
Using now pT|) . we deduce that also holds true. 



The barrier U^'^ is a subsolution of ^ on [to, to + h] xR^. We distinguish 
three cases. 

Consider first a point (t, x) such that x) < 7/2. In this case, U^'^{s, y) = 
U{s, y) in a neighbourhood of (i, x) and this implies dtU^'^{t, x) = dtU{t, x) (in 
the viscosity sense). In order to prove that W^'^ is a subsolution of (O at {t, x) it 
is enough to prove that l^U{t, x) < l^U'^-^{t, x) since U is a subsolution. Such 
an inequality is a consequence of ([60)1 and the fact that U'^'^{t, x) — U{t, x). 

Consider next a point {t,x) such that d{t,x) > 87/4. In this case, there 
exists ro > such that d{s,y) > 87/4 for y £ B{{t,x),ro). Consequently, 
W'l^is^y) = m+ for (s,y) e B{{t,x),ro). This yields that dtU"'/^ {t, x) = 
(in the viscosity sense) so we have to prove that 

f{U'^^'{t,x))<lUU'^^'){t,x). 

To get the previous inequality, on one hand, we have 

f{U'^''{t,x)) = /(m+-/3r/) 

< -Pv 
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and on the other hand, 

lUU'^^){t,x) = I [U'-'^{t,x + ez)- {m+~ P7j)]J{z)dz 

[U'''>^{t, x + ez)- (m+ - Pr])]J{z)dz 



\>ro/E 

> -C [ J{z)dz = -Ce" > -/3?7 

J\z\>ro/E 

in view of the definition of 77. Notice that this argument fails in the case a < 1. 

Finally, we consider {t,x) such that 7/2 < d{t,x) < 87/4. 

It is convenient to introduce ipdix) — ip{d{x)). Remark that U ~ U in a 
neighbourhood of {t,x). Hence, we mentioned above that U satisfies ([51]) at 
it,x) 

a,^ + l{_i^t/ + /(c7)}<^. (61) 

We use (pT|) and compute (in the viscosity sense) 



9tC/^"3 + — (/(C/^''')-X^[/^"3) 
er] 

^^'x (dtd) x{U~m+ + (ii^) + ijddtU + —fiW'l^) ~ —I%{U'-'^) 

erj erj 



U - {m+{-(3ri) - (iri) 



1 

erj 



i'd^ (62) 



where C^ip) only depends on ip and 7. We now estimate each term of the right 
hand side of ([5^ . 

First, we derive directly from the equality U = U and the very definition of 
U the following lemma 



Lemma 7. We have U = m+ — 2/377 + o{l3rj). In particular, 



We now estimate the second term of the right hand side of 
Lemma 8. 



1 



(63) 



(64) 



Proof. From Lemma [71 we have 



U'-^ = Vd(m+-2/3?7 + o(/3r7)) + (l-7/.d)(™+-/3??) 
= m+ - /377 - ■(/^^(l + o(l))/377. 
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In particular, U^'" < m+. Hence, 

fiU^'^) = .nU''P)-f{m+) + f{m+) 

< -{f'{m+) + o{l)){l + iPd + oil))(3rj-(3r, 

< -(/'(m+) + l + o(l))/3r? 

Lemma [7] also implies 

< U'^^ ~ U ^ {1 - + o(l))/3r; . 

Hence, we obtain 

< -(l-7/irf)(/'(m+) + l + o(l))/3r; 
+Mf'irn+) + o(l))(l - ^m){1 + o{l))Pv 

< (-l + o(l))(l-^,)/3,y<_i(l-V.d)/3r;. 

□ 

We now turn to the third term of the right hand side of (|62l) whose estimate 
is more delicate. It is given by the following technical lemma. 

Lemma 9. For any 70 > 0, there exists a function 0(e") such that 

^alltJ - TIU'^P < 7o'7 + 0(£") ■ (65) 

Proof. Wc would like first to point out that we can forget the time variable in 
this proof since it plays no role. 

We first remark that for ro = 7/4, 

|z| < ^ ^ If-'^ix + ez) = ^d{x + ez){lJ{x + ez) - m+ + l3r]) + (m+ - firj) . 

Indeed, \d{x + ez) — d{x)\ < tq ~ -f/A and this implies d{x + ez) G (7/4,7). In 
particular d{x + ez) < 7 and ([63|) holds true. 

We next approximate the quantity we are estimating by truncating large z's. 
Precisely, using the previous remark and the fact that the mass of J outside 
i?ro/e is 0(e"), we write 

iJdIlU -IIW'^ = ^dJ'U~J'[MU-m++l3rj) + {m+-Pi^)]+0{e'^) 
= ^dJ%U - m+ + /377) - J'['^d{U - m+ + /St?)] + 0(e") 

where the operator is defined as follows 

J^(p[x)— I [ip{x + ez) — ip{x) — D'^{x) ■ ezlB{z)\J{z)dz 

and where 0(e") only depends on HQIjoo and m± (for e small enough). 
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We next use the following equality 

{(p{x + ez) — (p{x)){ip' (x + ez) — ip' {x))J{z)dz 



f e\z\<rQ 

with ip ~ ^pd and ip' = U ~ m_|_ + fSrj. We obtain 

^bdl^U - T^U''" = -{U - m+ + PT])J'^bd 

{^dix + ez) - iPd{x))iU{x + ez) - U{x))J{z)dz + 0(e") . 



I e\z\<ro 

Recalling that holds true for \z\ < e^^ro, we write 
iPdKU-T^U'''^ < 2Pr^\J'Mx)\ 

+ e\\DtlJd\\L'^{Bix,ro)) (^4DU\\L^^B(x,ro))^ \z\^J(z)dz 

+ e\\Dijd\\L^[B{x.ro))'^P^ / |z|J(z)dz + 0(e"). 

J E<\z\<rae-'^ 

We next estimate each term as follows. 

\J'Mx)\ < c(^d)£", 

\z\^J{z)dz < Ce^-", 

z\<£ 

\z\J{z)dz < i r'\]-np\ ^fry~^ - 7' 

The first estimate is easily obtained by adapting the arguments used above to 
estimate R^[Tq] and R'^[Tq]. Moreover, the constant C{ipd) only depends on 
llV'dlloo — 1 and \\D'^ipd\\L'=°{B(x.ro)}- This last quantity only depends on 7 and 
||£'rf|lL~(B(a;,ro))i 1 1 -D^c^l I L~ (B(x,ro)) • Hcncc, by using Lemmaini we have 

ijdl'aU - r^U'^" < Cif^e" + e^-" + 77) + 0(£") 

(we used that (3 < 1 for instance). We achieve the proof by choosing e small 
enough so that 

C(r/£" + e^-" + 77) < 70 . 

□ 

We now combine (f62|) . ([64|) and ([65]) to get 
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This is where it is convenient to choose ip such that tp{5'y /8+r) = 1 — -0(57/8 — r) 
since in this case, max{'>pd, 1 — tpd) ^1/2 and we obtain 

Choosing now 70 small enough, we finally get 

er] 8e 

It is now clear that for e small enough, W^'^ is a subsolution of ^ in [to, to + 
h] X R^. □ 

5 Proof of Lemma [4] 

This section is devoted to the study of the average of oscillating terms. Their 
behaviour as £ ^ was given by Lemma |4] whose proof was postponed. We 
first deal with the singular case with a > 1. We next prove the result in the 
regular case. We then state the equivalent lemma for the case a < 1 since ideas 
will be used in the case a = 1. We finally prove Lemma 3] in the case a = 1. 

Proof of Lemma^ in the singular case for a > 1. We first recall the definition 
of and W. For the sake of clarity, we do not write e and h variables of q since 
they play no role in the present argument. 



ae{e,t,x)= J q{^)ae{^,e,t,x)dC 

= ^JJ qiO^Q{£. + e ■ z + eW(t,x,z)) - + e ■ z)y{z)d^dz 

with ^ 

W(t,x,z) = -ir [d(t, X + ez) — d(t, x) — eDd(t, x) ■ z] . 

£^ 

We proceed in several steps. 

Step 1: reduction to the study of the singular integral around the 
origin for quadratic Ws. Let us choose such that we also have (see ([S])) 



J{z)dz — > as £ — > . 

\>re 



For instance we consider r^ = e ^ with f3 > 1/a. In view of Condition ([5]), we 
thus can assume from now on that 
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Since q is bounded, it is therefore enough to study the convergence of 
bs{e,t,x) 

For \z\ < rs, 

W{t,x, z) ^D^d{t,x)z ■ z as £ ^ 

as soon as one chooses re such that er^ — > 0. Hence we take /3 € {cx~^, 1). If B 
denotes ^D'^d{t,x), we have for \z\ < rs, 

\W{t,x,z)- Bz- z\ <6\z\'^. 

By using the monotonicity of q, we thus can reduce the study of to the study 
of 

Ce = J q{0^^^q{^ + e■z + sCz■z)-q{^ + e■z)^J{z)d^dz . 

Step 2: integrating by parts. By using a system of coordinates where 
zi = e - z and 2; = (^i, ^;'), we can decompose the matrix C as follows 



C 

Hence, we can write 

\z\<r 



Cl V 

V C 



[ QiO^i^ + zi+ etCz ■ z){Cz ■ z)J{zi, z')d^dzidz'dT 

Cre Jo 

= [ [ [ q{£.)q{£, + zi+eTCz- z){Cz- z)J{z-i,z')d^dz-idz'dT 

J J\z\<r^ Jo 

= 11 f mdAq{^ + z^+eTCz-z)-q{£,)\ 



\z\< 

Cz ■ Z 

— J{zi, z')d£^dzidz'dT 



1 + 2eT{ciZi + v- z') 
We now integrate by parts with respect to Zi 



= -j j^^^ j\{0[<l{^ + zi+eTCz-z)^ 

{ 1 + 2eZ^+ V ■ z') Ji^^^^')]d^<^^^d^'dr + (BTy, - (BTt 
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with 



J J\z'\<reJ0 I 

+ £r(7(±Vr-f - \z'\^, z') ■ {±^rl - \z'\^, z')) 

L 1 + 2eT(±ci •s/rl - |z'|^ + ?; • 2;') J 
We next integrate by parts with respect to ^. 

+{BT)X - {BTy_ 
= 11 (\{0dz,\q{^ + z^+eTCz-z)-q{0\ 

J J\z\<reJo I J 

^ 1 + 2.r(crlr + . • .0 1 1 + 2er[^cr.r + . • z') ''^'^ ' ^oj^^^^i^^^^r 

We finaUy integrate by parts in z\ and we get 

c, = d, + {BT)\ - {BTf_ + {BT)\ - {BT f_ (66) 

with 

de = - I I f + ^1 + ^rCz ■ z) - q{0] 

J J\z\<reJO I J 

{ 1 + 2.r(crlr +v.z'Ml + 2erZz,\ v ■ z') ^(^i^ 

J J\z'\<reJ0 I 

+ST{±^rl - \z'\\ z') ■ (±Vr? - \z'\^ z')) - q{^)^ 



1 + 2£t(±ci - \z'\'^ +vz' 
Cz ■ z 



l + 2eT(ci(-) 



^ J(- , z') } (± Vr-I - md^dz'dr . 
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We now study the limits of all terms in ([M)) . 

Step 3: study of boundary terms. We start with {BT)j_. 

dz' 



\{BT)l\ < 2{m+{h) - m^{h))\\q\U\g\\^ I 2r, 



e N+a 



and this goes to as e ^ 0. 

We now turn to (i?T)j_. It is convenient to introduce the function 

r(r, zi) = (1 + 2£r(cizi + v ■ • 

Since er^ ^ as e ^ 0, we deduce that for e small enough, we have 

|r(r,zi)| < 2 

\d,,TiT,zi)\ < 8\ci\e 

|5.,(r2)(r,2i)l < 32|ci|e 

K,,,r(r,2i)l < 64|ci|V. 

We next compute 

a,,((Cz.z)r(zi)j(zi,z')) 

= 2cizirj + {cizf + C'z' ■ z'){d,,T)J+ {cizl + C z' ■ z')r9^i J . 
Now, since J{z) — g{z)\z\^^^°' , we deduce that 

\d.J{z)\< ^ 



N+a+l 



We thus conclude that for \z'\ < and zi such that Izl = r^, we have 



{{Cz ■ z)rj)\ < C-^ + + C-J^ < Crj^-'^ 

J. J. J. 



and we get 

\d,,{{Cz-z)r{T, zi)J{zi,z'))\dz' < CrJ". 



\z'\<r^,z\ + \z'\'^=rl 

With this inequality in hand, we now derive 



< {BT)l < Cr-'^ j {q{^ + Cr,) - qiO)d^ = G{Cr,) 



where 

G(r)= I iq{^ + r) - qiO)dC ■ 
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It is clear that G is Lipschitz continuous and equals at 0. Hence 

\{BT)l\<Crl-'^. 

It thus goes to as e — > 0. 

Step 4: study of d^. In order to study the main term d^, we first write it as 
follows 

4 = + Re 



with 



Re = - 



I I fq{0\(l{^ + Zi+eTCz-z)-q{0] 

J J\z\<re JO I J 

xr2(T, ^1)92^^^^ |(C^ • z)J{zi,z')^d^dzidz'dT 
[ [ f\{0{<li^ + zi+eTCz-z)-q{0] 

J J\z\<re Jo I ) 

I {d,,Vf + Tdl^^s] [ i.Cz ■ z)J{zu z')^ 



+d,,{T^)d,A{Cz-z)J{z,,z') 



d£,dzidz'dT . 



Let us now prove that Re goes to as £ ^ 0. We proceed as we did with 
(ST)j_. We first estimate the quantity [. . .] in the definition of -Re. We use 
the estimates on F and its derivatives, together with the estimate of d^^J. We 
obtain for \z\ < r^, 



|[---]| < c- 



c- 



< 



Ce 



|_2-|JV+a-2 ' l^-jiV+a-l — |_2,|iV+a-l 

since £|2:| < eVe — > 0. By arguing as for (BT)j_, we conclude that 
\Re\ < CeG{Cre) [ \z\-^-''+'^dz < C(£re)re^-" 

J\z\<re 

and the right hand side of the previous inequality goes to as e — > 0. 

It remains to study the limit of e^. By dominated convergence theorem, we 



30 



obtain that it converges towards 
eo = - [ q{0 I m + zi)-q{e}) 



+ ^i) - (Tr(A,(e)C7) 

where Ag(e) is defined in (f24|) . The proof is now complete. □ 

Proof of Lemma^in the regular case. The proof of the lemma in this case is 
divided into two steps. 

Step 1: reduction to the study of the singular integral around the 
origin for quadratic Ws. As explained in the proof of Lemma [4] for J{z) = 
g{z)\z\^^^" and a > 1, it is enough to study the convergence of 

J J\z\<r^ 

qiO-lqi^ + e • z + eW{t, X, z), e)-q{^ + e- z) \j{z)d£_dz 

l\z\<r, £ I J 

where we recall that ~ e^^ with (3 > 1/a. 

Remark that there exists Cu > such that for any {t, x) G Bj^ and z E B, 

\W{t,x,z)\<Cw\z\^ (67) 

and that, for \z\ < r^, 

W{t,x,z) ~> ■^D'^d{t,x)z ■ z as £ ^ 

as soon as one chooses r^ such that er^ 0. We conclude that the integrand 
of converges towards 

y'>{Oq'>{( + e-zKD'd{t,x)z-z) J{z). 

This explains why we expect the limit of a^ to be given by (|26[) . 

We can apply dominated convergence theorem outside the unit ball B. 
Hence, we reduce the study of the limit of to the one of 



Ce{e,t,x) 



l\z\<l 



q{i)-{ q{S, + e ■ z + eW{t,x, z),e) ~ q{£, + e ■ z) '>J{z)d^dz 
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In order to do so, we introduce 

- ' I^^^J{0l{^q{^ + e-z±eCw\z\\e)-q{^ + e-z)^J{z)d^dz. 

We know from (|67p and the monotonicity property of q that 

c7 <Ce< c+ (68) 

and it is thus enough to prove that integrals have hmits that are uniform 
with respect to e, t, x to conclude. 

Step 2: integrating by parts. Recall that |e| — 1 and let zi denote e • z and 
z — (zi, z'). We now write 



c+(e) 



9(0-^ 9(C + -zi + eCwzl + eCw\z'\^ ,e) ~ + zi) \J{z)d£,dz 



z\<l 



1 

2 , 1^/12 „\ I |2 



= Cw I / / + zi +eTCwzf + eTCw|zT,e) |zrJ(z) d^dzdr 

J J|2|<1 Jo 

^ 1 , o ^ — d£,dTdz'dzi . 

We next integrate by parts and get 

4 =dt + {BT)+ + (BT)^. 

with 

1 

/|2 



J Jq J\z\<1 

xK(t,z) d^drdz'dzi 

{BT)± = ±Cw I f I m<l{^±^^-\A^+eTCw.e) 
J Jq "'|2|<i 

^ J{±^/T^]lf,z') 



l±2eTCw\/l - 



where 



K{t,z) = d, 



1 + 2£tCwZi J 



1 + 2eTCH/zi (1 + 2£tCwziY 
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After remarking that for J E C] C^, we have 

\Kir, z)\ < Ci\z\'Jiz) + |V(|z|V(z))|) G L\B) , 

it is clear that we can apply dominated convergence theorem in each integral. 
The proof is now complete. □ 

In the case a < 1, the ansatz used to treat the case a > 1 yield oscillating 
terms with the following form. 

N ^ [ ( (/- , 4'is,y + £z)-<j){s,y) \ dz 

ae(C,s,y) = [q{C + : )-9(C + e-z)' 



J V^'^ ' £ y ■.v^ — 

Their average is thus defined as follows 

aeis,y) = / a^{r,s,y)q{r)dr . 



Even if we are not able (yet!) to treat the case a < 1, we think this can be of 
interest to explain what is the limit of the average as £ —> in order to justify 
our conjecture about the limit we expect in the case a < 1. Another reason for 
including such a result is that its proof shares ideas with the one corresponding 
one for the case a — 1. 

Lemma 10 (Uniform convergence of approximate coefficients (II)). Consider 
a smooth function d : (0, +oo) x M. such that \Dd{t,x)\ = 1. Then, as 

£ ^ 0, 

aj(e,i,x) — > K[x,d{t, •)] 
and the limit is uniform with respect to {e,t,x) £ E>^^^ x Qj. 

Proof of Lemma \l(A In the case a < 1, we first make a change of variables as 
follows 

de{e,t,x)^ I q{^)a^{^,e,t,x)dC 



^ f f I d{t,x + ez)~ d{t,x) \(Tf\^c^ 
liOs m + >e) -q{^ + e- z) \J{z)d£,dz 



£ _ 

with 

Step 1: reduction to the study of the singular integral around the 
origin for quadratic W^'s. Remark next that it is easy to pass to the limit in 
the integrand; indeed, 

d{t,x + z) - d{t,x) e-z 
q[i H , e - (7(C + 

£ £ 
-> (to+ - m_) {'i-(d{t,x + -)>d(t,x),e-(-)<0) - '^{d{t,x+-)«l(t,x),e-(-)>0)) 
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and 



N+a ■ 

The difficulty is to deal with the singular measure. Hence, it is enough to study, 
as in the case a — \, 

z\<S I £ ^ 

X (zi , z')d£^dz'dzi 

JMf J\z\<S Jo ^ 

X Js{zi, z')d^dzdT 

JMf J\z\<S Jo \ ^ 

where z — z\e + z'. 

Step 2: integrating by parts. By integrating by parts, we obtain 

c+ = rf+ + (Br)+ + {BT)_ 

where 

jRj J|g|<(5 Jo £ 

xLg{T, z)d^dzdT 

inn, = ±c,,/ / /',K„«.M!ESi±^w!^^, 

JKf J|z'|<<5 Jo £ 



c<5Ve(±v/(52- z')dz'd^dT 



with 

1 + 2Cvi/TZi (1 + 2CwTZi)^ ■ 

Condition © ensures that 

\L,{T,z)\<CK{z)eL\B) 



Le{T,z) = 
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and dominated convergence can be used to prove the convergence of . As far 
as boundary terms are concerned, we use ([8]) to get a constant C > such that 



|7V+Q 

and this imphes 

Hence, dominated convergence can be appHed to boundary terms too. The proof 
is now complete. □ 

Proof of Lemma^in the singular case for a = 1. The proof is divided in sev- 
eral steps. 

Step 1: reduction to the study of the singular integral around the 
origin for quadratic Ws. Let us fix 5 > 0. There exists rs > such that for 
any z e Br^, 

\d{t, x + z)- d{t, x) - Dd{t, x) ■ z - ^D'^d{t, x)z ■ z\ < S\z\'^ . 
Consequently, for z such that \z\ < rse^^ i?^, we have 

\W{t,x,z)-^Bz- z\<S\zf . (69) 

We remark next that 

\z\^J{z)dz = 0(1) = o(|ln£|), 

z\<l 



[ J{z)dz = 0{e) = o{e\\ne\ 



It is therefore enough to study the convergence of 



e\\ne\ J Ji<\,^\<r 



q{Cj[(l{S. +e- z + eW{t, X, z)) - + e ■ z)]J{z)dzd£, 



By using together with the monotonicity of q, it is even enough to study 
the convergence of 



1 



e|lne| JJi<\z\<R. 



QiOki^. +Z1+ ecizj + eC'z' ■ z') - + e ■ z)]J{z)dzd^ 



for any {N— 1) x {N— 1) symmetric matrix C and any constant ci £ K. Precisely, 
we next prove that, as e ^ 0, 

c, 2(m+ - m_)^ / g{e) C'O ■ 9 a{de) 
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where g appears in ([5]). 

Step 2: change of variables and domain decomposition. We now intro- 
duce the function Fq defined as follows 



We remark that Fq is bounded, non-decreasing, Lipschitz continuous and sast- 
isfies 

m+(m-|- — m_) as a ^ +00 
m-lfh+ — ffi-) as a ^ —00 



and 



1/2 



We now rewrite Ce with this new function 
1 



el Inel 



Fq{zi + eciz{ + eC'z' ■ z') - Fq{zi) \J{z)dz 



Through the change of variables z' = e ^rO and zi = e ^rti, we get 

1 r dr r 

da{0) / ^ / dh 



line 



Fq (^(ti + rcitj + rC'6 ■ 0)] -Fq (^ti) 



Without loss of generality, we assumed here for simplicity that 

Jiz)=g{z/\z\)\z\-^-' 
(see ©)■ We next cut it into two pieces as follows 

Ce = + Rl 

with 

1 /■ /"''' dr f 

de = ^ r / da{e) — dti 

|lne|JgggN-2 f Je^r-^<l+tj<r]r-2 

x\Fq{...)-Fq{...)}g{...)- ^ 

Rl - 



(l+t2)(JV+l)/2 ■ 

— f dai0)r^f dt. 



y.\Fq{...)-Fq{...)\g{...) ^ 



(l+t2)(JV+l)/2 
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for some 6 > to be fixed later. 

Step 3: study of i?^. Wc now use the fact that Fg is Lispchitz continuous in 

order to get, for 6 large enough, 

1 r^^ 

\Rl\ < IMIcx>||F^||oo(|ci| + |C"|)|S^-2|^y^ dr 



L 



, 1 1 

dti-- 



5 



£| inS^I Jo Jl<|ti|<rir-1 



Hence 



Step 4: study of de- First, we rewrite d^ (for 6 > 1 large enough) as follows 



I m£| Jq 



1 



X / , {F,i...)-Fgi...)]g{...) 



We now study the limit of the integrand of d^- In view of the limits of Fg{a) 
when a ±oo, this limit is not if ti satisfies 

ti>0 ( ti<0 

h + cirtl + rc0<O °^ \ h+ cirtl + rcg > 

where ce = C'6 ■ 9. This is equivalent to 

0<t.<T,(r) or ( < "dr^d + VT^^^) - 
[or — cer < < 

with ^ 

TlM = - — - 4cicer2 - 1) ~ -C6ir . 
zcir 

We use now the fact that < -\/r|r~2 — 1 to get 

< ii < Ti(r) or Ti(r) < ii < . 

We next cut d^ into pieces as follows 

4 = e+ + er + 
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with 



4 = Ti^l da{e)r% 



I In 

x\F,{...)-F,{...)}g{...) 



eeB"-^ Jes ' Jo<±ti<\Ti{r)\ 

1 



(l+i2)(JV+l)/2 



(■r)\<\tl\<^/rjr-^-l 
1 

'(l+i2)(JV+l)/2 ■ 

Step 5: study of e^. Since r < rs and very small, we can replace ef with 
I inel J9gsN-2 r Jo<±ti<|ri(r) 

dr f 

da (9) — dzi 

eGS"-2 Joe JO<±zi<r-i|Ti(r)| 



llnel 



Now it is easy to conclude that 



lee 



I ((C0)+ + {ce)-)g{Q,e)da{e){F^o{+^) - F,p{~^)) 



(m+ - m-Y / (Ce • %(0, 6')dCT(6') 



It remains to prove that converges towards 0. 

Step 6: study of R^. By the study we made in Step 4, we conclude that the 
integrand of R^ converges towards 0. To dominate convergence, we simply write 

\Rl\ < T- / 2||F, 



Inel' '7ee r 7|.,|>|c.|/2 " -^""^ (1 + tf 
< |§A^-2| ln(^^) + |lne| f ^HF,!' ^ 



|lne| 7|.,|>|.«|/2 " ^"°°(l+t2)(^+i)/2- 
The proof of the lemma is now complete. □ 



A Link between O and (I2SD 

In this section, we explain the link between the two first formulae in (|23p 
and (|26p by giving a formal argument. 
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If zi denotes e • z and z' = z — [e ■ z)e G ^, then this equation can be 
written as follows 



Aie) = 



[ d^q'iO [ dziq%^ + zi) 
We remark next that 

Zf Z\Z 

zi{z')* z'(g)z' 
Hence, A{e) = KAg{e) with 

K ' 



dz' 



J{zi,z')dz' 



Zl Z\Z 

zx(z')* z'®z' 



Ag{e)\zx 



J{zi,z') 



\l-a 



a{a — 1) 



, d^q"{0 / dz,qO{^ + z,)\z,\'-". 



(70) 



By integrating by parts in zi and ^ and zi successively, we obtain 

Ji = - [ d^qO{0 [ dz,{q°{^ + z,)-q°mzir^'^''^zi 

= -- / diq\i)\ d-2i(g0(^ + ^i)-g0(0)kir<'+"^^i 
= -[ d^q°{0 [ dz,{q°{^ + z,)-q°{0-q°{Ozi)\zi\-'' 

= J g°(0(-A)(«-^)/^[g°](Orf^ 

= / J{A^ + zi)-Aorj^d^- 
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